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1. Introduction 

Much work has been devoted in the last years to the study of the QCD running coupling 
constant determined from lattice simulations, as well in its perturbative regime [1-9] as in the deep 
infrared domain [10]. Only very recently [9, 1 1], the Green's function approach to study the running 
coupling and then to estimate A^jg- has been pursued by exploiting a non-perturbative definition of 
the coupling derived from the ghost-gluon vertex. The very infrared domain for the running of the 
coupling so defined has been discussed in the Olivier Pene's talk. We aim to deal here with the 
running behaviour of this ghost-gluon coupling beyond the IR domain, above roughly 2-3 Gev. 

We will show that the analysis of quenched lattice simulations leads to a non-perturbative 
determination of the running coupling in terms of two-point ghost and gluon Green functions and 
to obtain A^r in pure Yang-Mills (Nf = 0). Furthermore, a very realistic estimate of A^, directly 
comparable with experimental determinations, will become an immediate possibility thanks to the 
many unquenched configurations which are presently available. 

A precise determination of the non-perturbative coupling from the lattice also reveals a dimension- 
two non-zero gluon condensate in the landau gauge [10]. One needs then to describe the running 
with a formula including non-perturbative power corrections to be confronted with lattice estimates 
of the coupling. This procedure constitutes an optimal method for the identification of A^ and of 
the gluon condensate [11]. Much work has been also done to investigate its phenomenological im- 
plications in the gauge-invariant world [12]. In particular, we will discuss the interpretation of this 
condensate in terms of the Yang-Mills semiclassical field background by applying the Instanton 
liquid model. 

2. The ghost-gluon coupling 

There is a large number of possibilities to define the QCD renormalized coupling constant, 
depending on the observable used to measure it and on the renormalization scheme. Actually, any 
observable which behaves, from the perturbative point of view, as g provides a suitable definition 
for it. Among such quantities stand the 3-gluon and the ghost-gluon vertices, which have been 
widely used by the lattice community to get a direct knowledge of a s from simulations. Of course 
an important criterion to choose among those definitions will be how easy it is to connect it to other 
commonly used definitions, specially the MS one, and to extract from it fundamental parameters 
like A QCD . 

A convenient class of renormalization schemes to work with on the lattice is made of the so- 
called MOM schemes which are defined through the requirement that a given scalar coefficient 
function of the Green's function under consideration take its tree -level value in a specific kinemat- 
ical situation given up to an overall "renormalization scale". To make the point clearer we recall 2 
schemes which we have used in previous works on a s : 

• The symmetric 3-gluon scheme in which one uses the 3-gluon vertex T llV p(pi,p2,P3) with 

2 2 2 2 

• The asymmetric 3-gluon scheme (MOM) in which the 3-gluon vertex T^ vp (pi,p2,P3) is 
used with p\ = p\ = jx 2 , p\ = 



2 
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In the present note we shall apply a specific MOM-type renormalization scheme defined by 
fixing the (ghost and gluon) propagators and the ghost-gluon vertex at the renormalization point. 
Let us start by writing the ghost and gluon propagators in Landau gauge as follows, 



P>)V,A) = ^^(V£^) 
/fiv p 1 \ P J 



^ 2a , s F(p 2 ,A) 



' (/,A) = -S a »^^; (2.1) 



P 2 

A being some regularisation parameter if, for instance, we specialise to lattice regularisa- 

tion). The renormalized dressing functions, Gr and Fr are defined through : 

G R (p 2 ,^ 2 ) = lim A ^Z 3 - 1 (Ai 2 ,A)G(^ 2 ,A) 

F R (p 2 ,^ 2 ) = \im^ 00 Z-\pL 2 ) K)F(p 2 A) , (2.2) 
with renormalization condition 

G s (At 2 ,At 2 )=^(M 2 ,At 2 ) = l • (2-3) 

Now, we will consider the ghost-gluon vertex which could be non-perturbatively obtained through 
a three-point Green function, defined by two ghost and one gluon fields, with amputated legs after 
dividing by two ghost and one gluon propagators. This vertex can be written quite generally as: 

ff c (-q,k;q-k)= - ^ " = igof^ (q v Ht(q,k) + (q - k) v H 2 (q,k)) , (2.4) 
k q 

where q is the outgoing ghost momentum and k the incoming one, and renormalized according to: 

f R = Z l r. (2.5) 

The vertex T v involves two independent scalar functions. In the MOM renormalization procedure 
Z\ is fully determined by demanding that one specific combination of those two form factors (cho- 
sen at one's will) be equal to its tree-level value for a specific kinematical configuration. We choose 



to apply MOM prescription for the scalar function H[ + H2 that multiplies q v in eq. ( |2.4[ ) and the 
renormalization condition reads 1 

(H?(q,k)+H«(q,k))\ q2= ^ 2 = Urn Z^ 2 ,A) (H^k; A) +H 2 (q,k; A)) ^ = 1, (2-6) 

where we prescribe a kinematics for the subtraction point such that the outgoing ghost momentum 
is evaluated at the renormalization scale, while the incoming one, k, depends on the choice of 
several possible configurations; for instance: k 2 = (q — k) 2 = \i 2 (symmetric configuration) or 
k = 0, (q — k) 2 = }X 2 (asymmetric-ghost configuration). 



'in the case of zero-momentum gluon, an appropriate choice would be Z\ ([l 2 )H\ (q.q)]^^^! = 1. This would make 
the renormalized vertex equal to its tree-level value at the renormalization scale. 
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On the other hand, the fields involved in the non-perturbative definition of the vertex T v in 



eq. ( |2.4| ) can be directly renormalized by their renormalization constants, Z3 and Z3, and the same 
MOM prescription applied to the scalar combination H\ +H2 also implies: 



lim Z3(Ai 2 ,A)Z 3 1/2 (At 2 ,A)^o(A 2 ) ( #iM;A) +H 2 (q,k;A) 



Zi(M 2 ,A 2 ) 



lim g (A 2 ) 



(2.7) 



We combine both eq. ( |2.6| ) and the first-line equation of ( |2.7| ) to replace H\ + H 2 and obtain the 
second line that shows the well-known relationship Z g = (z\ Zt,)~ 1 Z\, where gn = Z~ l gQ. 

We turn now to the specific MOM-type renormalization scheme defined by a zero incoming 
ghost momentum. Since those kinematics are the ones (and the only ones) in which Taylor's well 
known non-renormalization theorem (cf. ref [13]) is valid we shall refer to this scheme as to the 



T-scheme and the corresponding quantities will bear a T subscript. Then, in eq (p.4|), we set k to 
and get 



ff c (- 



q,0;q) =igof bc (H l (q,0)+H 2 (q,0)) q v 



(2.8) 



Now, Taylor's theorem states that H\(q,0;A) +H 2 (q,0;A) is equal to 1 in full QCD for any value 
of q. Therefore, the renormalization condition eq. ( |2.6| ) implies Zi(/i 2 ) = 1 and then 

S 2 T (H 2 ) 



OrOi ) = 



An 



lim 



gg(A 2 ) 
An 



G(ii 2 ,A z )F z (Li z ,A' 



(2.9) 



where we also apply the renormalization condition for the propagators, eqs. ( |2.2| , [2.3| ), to replace 
the renormalization constants, Z3 and Z3, by the bare dressing functions. The remarkable feature 



of eq. ( |2.9j ) is that it involves only F and G so that no measure of the ghost-gluon vertex is needed 
for the determination of the coupling constant. 



Equation (2.9) has extensively been advocated and studied on the lattice (see for instance 
reference [14]). However it must be stressed that the T-scheme is the only one in which Z\ = 



1. Nevertheless the form (2.9) is used quite often in this case (for a kinematical configuration 



other than T-scheme's) also as an approximation, specially in relation with the study of Dyson- 
Schwinger equations. An important remark is also in order here: in the very infrared domain, for 
phenomenological purposes (see for instance [15]), the coupling can be more properly defined by 
pulling a massive gluon propagator out from the ghost-gluon Green function used to build it [16]. 

2.1 Pure perturbation theory 

A standard four-loop formula describing the running for the T-scheme coupling, 



orQfi) = £ , ft !2l« + g _iy + « ' 



+ 



Pot 



1 



PI 

ft 



+ ■ 



1 



(ftf) 4 \2Po 2 VA) 
.,2 




P\ 



2 log 3 (0 + 5 log 2 (0+ 4-6 



P2P0 
Pf 



logO) - 1 



with t = In -^r 
Aj 



(2.10) 
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is obtained by inverting the f5 -function of a,j, 

where, as explained in [11], the coefficients j3 ( can be computed in terms of j3 ( -, those for the j3- 
function of the coupling renormalizad according MS-scheme, a, and of the anomalous dimensions 
for gluon and ghost propagators, 

1 da T (u 2 ) 1 /„ ,. d , „. o A . d , „. o A . 

'2hm— T lnF jU 2 ,A + hm — T lnG(/x ,A) 



ar(;U 2 ) <ia AtfsC 01 ) V A ^°°dlnn 2 ' A^~<iln/i 

2y(a) + y(a 
~ /W«) 



(2.12) 



Both anomalous dimensions need to be renormalized along MOM prescriptions {i.e., Gr(/j. 2 ,/j. 2 ) = 
Fr(/j. 2 , n 2 ) = 1) but expanded in terms of a. The coefficients so obtained (the details and results 
of the computation can be found in [11]) appear to agree with those directly obtained in ref. [17] 
by the three-loop perturbative substracion of the ghost-gluon-gluon vertex in the QCD Lagrangian 
with the appropriate kinematical configuration (J-scheme). 

2.2 OPE power corrections 

In order to extend the description of the running coupling down to energies as low as possible 
(of the order of 3 GeV) and to take full advantage of the lattice data we want to compare with, 
in order to reduce the systematic uncertainties, it is mandatory to take into account the gauge- 
dependent dimension-two OPE power corrections (cf. [7, 8, 10, 18]) to aj- 

The leading power contribution to the ghost propagator, 

{F [2) ) a \q 2 ) = J d 4 xe iqx ( T (>(*)?(())) ) (2.13) 
can be computed using the operator product expansion [19] (OPE), as is done in ref. [20], 

T (c a (xy(0)) =£(c t ) ai (x) O,(0); (2.14) 

here O t is a local operator, regular when x — > 0, and the Wilson coefficient c t contains the short- 
distance singularity. Eq. ( [2.14 ) involves a full hierarchy of terms, ordered according to their mass- 



dimension, among which only 1 and : A a ^A v : contribute to eq. ( [2.1 3\ > in Landau gauge up to the 
order l/q 4 . Then, using eq. (2.14) into eq. (2.13), we obtain: 



(F^V) = (c rV) + (c 2 r s ^(q 2 )(:AU0)K(0)--) + ••• 

= (4*rV) + ™ ab 4{ ^ 2) _ T) + ■■ (2.i5) 



2 The operators with an odd number of fields (d = 1,3/2; d^A and d^c) cannot satisfy colour and Lorentz invariance 
and do not contribute a non-zero non-perturbative expectation value, and cc does not contribute either because of the 
particular tensorial structure of the ghost-gluon vertex. 
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where 



w 



,ab 



2x 




Jd 4 xe^ x (A f T ',(0) T (c a c h ^ A^(0)) 



connected 



(2.16) 



and the SVZ factorisation [21] is invoked to compute the Wilson coefficients. Thus, one should 



compute the "sunset" diagram of the last line of eq. ( J2. 16| ), that binds the ghost propagator to the 
gluon condensate (where the blue bubble means contracting the color and lorentz indices of the 
incoming legs with \/28 st 8 ax ) to obtain the leading non-perturbative contribution (of course, the 
first Wilson coefficient gives trivially the perturbative propagator). 

Finally, after the renormalization of the A 2 -condensate at the subtraction point q 
cording to the MOM scheme definition, the ghost dressing function is written as: 



H 2 , ac- 




(2.17) 



where the multiplicative correction to the purely perturbative Fs.pert is determined up to corrections 
of the order l/q 4 or Inq/jj,. As far as we do not deal with the anomalous dimension of the A 2 



operator, factorising this purely perturbative ghost dressing function in eq. ( |2.17[ ) is a matter of 
choice. However, the Wilson coefficient is also computed at the leading logarithm in ref. [11] and 



eq. (117) appears then to be a very good approximation up to this order. 

We can handle in the same way (see refs. [7,8]) the OPE power correction to the gluon propa- 
gator, 



pert/JivV 



z 2 ,M 2 ) 



+ [W 



ab 
" [IV 



(A 2 



Ry 4(a£- l) 



+ 



(2.18) 



and obtain 



w 



.ab 




+ 2 




3 £ 2 (r (2) Y b 

2 V-^perU/iv • 

H 



(2.19) 



Then, after renormalization and appropriate projection, one gets for the gluon dressing function: 

J2//t2\ 



G R (q 2 ,H 2 ) = G R 



pert \q i 



4 ' } { Y + q 2 4(N 2 -l) J • 



(2.20) 
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Finally, putting together the defining relation eq. (2.9) and the results eqs. (2.17, 2720| ) we get 



arGu 2 ) = lim ^-F 2 (p. 2 ,A)G(p- 2 ,A) 



lim ^-F 2 (q 2 ,A)G(qt h A) F^(^,q 2 ) G R (^,q 2 ) 

,pert/2 W 2 ,,.2 „2^ ^ ( ..I JL\ ( <i , 9 (A 2 ) S , 9 2 



of (?o)^ P ertG" 2 ,<7o) GrMV i<ll) 1 + 



M 4(a£-i; 



(2.21) 



where g 2 , 3> Aqcd is some perturbative scale and the running of the perturbative part of the evo- 
lution, a^ ert , is of course described by the eq. ( 2. 1C ) in the previous section. Again, the Wilson 
coefficient at leading logarithm for the T-scheme MOM running coupling is obtained in [11] and 
found not to induce a significant effect, provided that the coupling multiplying A 2 inside the bracket 



is taken to be renormalized also in T-scheme. Thus, eq. (2.21) describes pretty well the running of 
(Xt roughly above 3 Gev. 



3. Data Analysis by the "plateau" method 



In the following, as done in [11], we will apply a "^/afeaw"-procedure exploiting eq. ( [2.2 1[ ) to 
get a reliable estimate of the AQCD-parameter from lattice data. The goal being to get a trustworthy 
estimate of the A^-parameter, one could attempt to do it by inverting the perturbative formula 



eq. (2.1C) and using in the inverted formula the lattice estimates of the running coupling obtained by 



means of eq. ( |2.9[ ) for as many lattice momenta as possible. Then, one should look for a "plateau" of 
Aj^jg in terms of momenta in the high-energy perturbative regime (this was done with the coupling 
defined by the three-gluon vertex in [4,5]). In the next subsection, fig. |j].(a) shows the estimates of 
T^ms so calculated for the lattice data presented in ref. [20, 22] over 9 <p 2 < 33 GeV 2 . 

However, in order to take advantage of the largest possible momenta window one can use 



instead eq. (|2.21| ). In this way we shall hopefully be able to extend towards low momenta the 



region over which to look for the best possible values of the gluon condensate and of Aj^r . In 
other words, one requires the best-fit to a constant of 



(pf,A(ai)) , with: a,- 



«Latt(^ 2 

i+4 



(3.1) 



where A(a) is obtained by inverting the perturbative four-loop formula, eq. ( |2. 10h , and c results 



from the best-fit (it appeared written in terms of the gluon condensate in eq. ( |2.21[ ) ). Thus, A(a) 
reaches a "plateau" (if it does) behaving in terms of the momentum as a constant that we will take 
as our estimate of A^. Of course, this is nothing but a fitting strategy for a 2-parameters (A^ an d 
(A )) fit of the estimates of eq. ( ^9| ) from lattice data. 



3 This increases the statistics and reduces errors. It also avoids some possible systematic deviation appearing when 
lattice momentum components, in lattice units, approach n/2 (Brillouin's region border). 
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3.1 Results for pure Yang-Mills (N f = 0) 

The quenched lattice data that we will exploit now were presented for the first time in ref. [22]. 
We refer to this work for all the details concerning the lattice implementation: algorithms, action, 
Faddeev-Popov operator inversion, etc. The parameters of the whole set of simulations are de- 
scribed in table [I] 



p 


Volume 


a' 1 (GeV) 


Number of confs. 


6.0 


16 4 


1.96 


1000 


6.0 


24 4 


1.96 


500 


6.2 


24 4 


2.75 


500 


6.4 


32 4 


3.66 


250 



Table 1: Run parameters of the exploited data. 



In fig. 0.(a), we show the estimates of A^ obtained when interpreting the lattice coupling 



computed by eq. (2.9) for any momentum 9 <p 2 <33 GeV 2 in terms of the inverted four-loop 



perturbative formula for the coupling, eq. ( |2. 10[ ). The estimates systematically decrease as the 
squared momentum increases until around 22 GeV 2 ; above this value, only a noisy pattern results. 
In fig. |l|.(b), the same is plotted but inverting instead the non-perturbative formula including power 



corrections, eq. (2.21). The value of the gluon condensate has been determined by requiring a 
"plateau" to exist (as explained in the previous section) over the total momenta window. 



0.28 



0.22 





(a) 



(b) 



Figure 1: (a) Plot of A-^g (in GeV) computed by the inversion of the four-loop perturbative formula as a 
function of the square of the momentum (in GeV 2 ); the coupling is estimated from the lattice data through 
the perturbative formula obtained in the text, (b) Same as plot (a) except for applying the non-perturbative 
formula obtained in the text for the coupling and looking for the gluon condensate generating the best plateau 
over 9 < p 2 < 33 GeV 2 . 

One should realize that the non-perturbative analysis seems to indicate that the perturbative 
regime is far from being achieved at p = 5 GeV. This is also illustrated by figure |2|.a in which, 
adopting for A^- the value 224 MeV which results from the non-perturbative analysis, we plot 
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F 2 G [11] 


Asym. 3-g [8] 


Sym. 3-g [8] 


F/G [20] 


[2] 


AMs(MeV) 


224+* 


260(18) 


233(28) 


270(30) 


238(19) 


V(A 2 ) RjP (GeV) 


1.64(17) 


2.3(6) 


1.9(3) 


1.3(4) 





Table 2: Comparison of the estimate of A-^g obtained from the analysis of the ghost-gluon vertex (first 
column) and others from literature. The renormalization point is fi = 10 GeV. 



against the square of the renormalization momentum the coupling constant as computed by means 
of the non-perturbative formula (2.21) (red curve) and of the perturbative one ( 2.10| ) (blue curve). 
Displayed are also the lattice data, i.e. the values of otj obtained from eq. (E 



Average 




(a) 



(b) 



Figure 2: (a) Plot of <Xj in terms of the square of the renormalization momentum: the red solid line is 
computed with the non-perturbative formula with Aj^g = 224 MeV, the blue one with the perturbative one 
for the same A^g and the data are obtained from the lattice data set-up described in the text, (b) Comparison 
with previous published estimates of A^g in pure Yang-Mills; the blue triangle stands for the estimate in this 
work and the red square for the average of the five estimates presented in the plot. 

Thus, one can conclude that our best-fit parameters incorporating only 4 statistical errors are: 

A^T° = 224+ 8 MeV 

MS 

g 2 T {A 2 ) R = 5.1+^GeV 2 . (3.2) 

These values are in very good agreement with the previous estimates from quenched lattice simu- 
lations of the three-gluon Green function [7, 8] or, in the case of A^g, from the implementation of 
the Schrodinger functional method [2], although slightly larger than the one obtained by the ratio 
of ghost and gluon dressing functions [20] (see fig. |2[(b) and tab. f|). 

4. About the nature and the size of the gluon condensate 

The nature of the dimension-two gluon condensate, as well as its possible phenomenological 
implications, have been discussed in many works in the last few years (see for instance [12, 18]). 

4 The error analysis is deeply discussed in [11]. 
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In particular, we presented some indications supporting the idea that the low-momentum gluon 
correlation functions could be nicely described in terms of the semiclassical instanton background 
for the gauge field [23], and used an instanton liquid picture to estimate the size for this gluon 
condensate in Yang-Mills [10]. Indeed, the gauge field in the instanton picture and within the 
sum-ansatz approach, can be written in the singular Landau gauge as 

=2l«v^^0 (^) , (4-1) 

where g = (6//3) 1 / 2 is the bare gauge coupling in terms of the lattice parameter j3, tj is known as 
't Hooft symbol and R aa represents the color rotations embedding the canonical SU(2) instanton 
solution in the SU(3) gauge group, a = 1, • • • ,3 {a = 1, • ■ • ,8) being an SU(2) (SU(3)) color index. 
The sum is extended over all the instantons and anti-instantons (we should then replace the 't Hooft 
symbol tj by tj) in the classical background of the gauge configuration. <p (x) is the instanton profile 
function. If we consider the profile of the BPST solution for an isolated instanton, we get 

g 2 <A 2 > = *±?± f ^£,A« 8 A« = 12*V^ = nnVn ; (4.2) 
where Nj (Na) stands for the total number of instantons (anti-instantons). On the other hand, if we 



neglect instanton position and color correlations, eq. ( |4.1| ) leads for the m-gluon Green function to 



Ak 2 ( R \ m l 2 

G< "°<* 2 > = (4) <p 3 w> 



where I(s) = / z,dz,J 2 {sz) <p(z) , (4.3) 

s Jo 

for m = 2,3; n being the instanton density. It depends on the functional I(kp) of the general 
instanton profile, (x), and < • • • > means the average over instanton sizes with a given normalised 
instanton radius distribution, ju(p). Then, two interesting limits appear where some results not 
depending on the instanton profile can be obtained: 

• For a sharp radius distribution, the particular combination of two and three-gluon Green 
functions defining the three-gluon running coupling in ref. [5] gives [24] 

k 6 (G^) 2 k 4 
a (jk 2) = ^^-J- = J— ; (4.4) 
gv ; 47T (g( 2 )) 3 18tt« 

• For kp S> 1, as I(s) asymptotically behaves as \6n 2 /s 3 in the large s limit, one obtains 

~ „i (^lY' 2 ^ (4 5) 

m \ 3 / 



Thus, eq. (4.5) provides us with large-momentum limits for the two and three-gluon Green func- 
tions behaviour which do not depend on the radius distribution nor on the instanton profile. How- 
ever, the large-momentum lattice correlation function being dominated by the short-distance quan- 
tum fluctuations, whether such a behaviour occurs can be only detected after performing some 
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"cooling" procedure [25] to kill the higher energy modes. This is done in [23] and, as can be seen 
in fig. |[a, the expected k~ 6 (k~ w ) power behaviour clearly emerges for the two-gluon (three-gluon) 
Green function after "cooling". This is a good indication for the success of the instanton picture 
in describing the gluon correlation functions. However, as the "cooling" has been proved to alter 
the configuration (instanton sizes become distorted, instanton and anti-instanton anhiliate to each 



other...), the power-law given by eq. (|4.4|), which is thought to be followed by the "uncooled" gluon 
correlators in the low-momentum regime, offers a more reliable "instanton detector". In ref. [24], 
eq. ( [O] ) is shown to work for a three-gluon coupling computed from several lattice simulations 
(see fig. ||.b taken from [24]) and the instanton density is estimated to be n ~ 5 fm -4 . 





0.25 0.5 
k (GeV) 



(a) 



(b) 



Figure 3: (a) two and three-gluon Green functions after cooling: they reach their expected power-law when 
the number of cooling sweeps increases, (b) The three-gluon coupling defined in the text: it follows the 
expected low-momentum k 4 power-law with n = 5.27(4) fm -4 . 

Then, this estimate of the instanton density and the average instanton radius, p ~ 0.4 fm 
(measured, for instance, in [23] and being close to the phenomenological prediction, ~ 1/3 fm), 



can be applied to eq. (<L2) to give: g 2 (A 2 ) ~ 4GeV 2 . There is of course no exact recipe to compare 
this estimate with the OPE one, since the separation between the semiclassical non perturbative 
domain and the perturbative one cannot be exact 5 . However, both lie prettily on the same ballpark. 



5. Conclusion 



We have demonstrated that, in the particular T-scheme, the coupling defined from the ghost- 
gluon vertex is obtained by only dealing with two-point Green functions. Some interesting non- 
perturbative information can be furthermore distilled from the running analysis of this coupling 
beyond the deep IR: the Aqcd parameter (usually expressed in the MS-scheme), computed here 
for pure Yang-Mills from Landau gauge lattice simulations, and a gauge-dependent dimension-two 
gluon condensate. The latter is interpreted and sized by invoking an instanton liquid picture, which 
successfully describes the low-momentum gluon correlations. 

5 One may appeal to the fact that at the renormalisation point fi, the radiative corrections are minimised; therefore 
a semiclassical estimate must best correspond to (A 2 )r ^ at some reasonable fi, which one could guess to be a typical 
scale of the problem as 1/p or some gluon mass. 
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